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Abstract. - We examine the effect of Coulomb interaction on the mobility and diffusion of 
spin packets in doped semiconductors. We find that the diffusion constant is reduced, relative 
to its non-interacting value, by the combined effect of Coulomb-enhanced spin susceptibility 
and spin Coulomb drag. In ferromagnetic semiconductors, the spin diffusion constant vanishes 
at the ferromagnetic transition temperature. 



The ability to control inhomogeneous distributions of electrons and holes in semiconductors 
is essential to the operation of modern electronic devices. Unlike electron packets in metals, 
which spread out very quickly under the action of their own electric field, electron-hole pack- 
ets in semiconductors are charge-neutral objects and can therefore be long-lived. The time 
evolution of such packets is controlled by a drift-diffusion equation with mobility and diffusion 
constants /j and D, as was verified in detail in the classic Haynes-Shockley experiment Q. 

Recently, a broader category of possible disturbances, involving inhomogeneous spin dis- 
tributions, has come into sharp focus in the context of the emerging field of "spintronics" S. 
Consider, for example, a spin packet consisting of excess up-spin electrons compensated by 
a deficiency of down-spin electrons. Such a disturbance can occur in the conduction band 
of a metal or of a doped semiconductor Like an ordinary electron-hole packet this is a 
charge-neutral object and can therefore be extremely long-lived (recent experiments measure 
a spin relaxation time t s of the order of 10 ns HQ). Unlike an electron-hole packet, however, 
the disturbance involves carriers of a single polarity - electrons - and therefore evolves with 
the mobility and diffusion constants of the conduction band, which are usually larger than 
those of the valence band. Indeed, large values of the spin drift mobility /i s ~3x 10 3 cm 2 /Vs 
and spin diffusion constant D s ~ 10 3 cm 2 / s have been recently observed in experiments on 
n-doped GaAs [Q. 

The qualitative difference between unipolar and bipolar disturbances in semiconductors 
has recently been emphasized by Flatte and Byers j^] within the frame of a simple model 
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in which the electron-electron interaction is treated in the Hartree approximation. In this 
Letter, we refine their analysis by examining the more subtle effects of exchange and Coulomb 
correlation on the mobility and diffusion constant of a spin packet. It will be shown that these 
many-body effects have a serious impact on the value of the spin diffusion constant and can 
be quantitatively probed in a Haynes-Shockley-type experiment that measures independently 
the mobility and the diffusion constant of a spin packet. The two key physical effects are (i) 
the reduction of the spin stiffness (inverse of the spin susceptibility) due to (mainly) exchange 
interactions and (ii) the spin Coulomb drag Q working like friction against the relative motion 
of up-spin and down-spin electrons. Both effects tend to reduce the diffusion constant. By 
contrast, the spin-packet mobility turns out to be essentially unaffected by interactions. 

Under certain conditions, the electron (hole) gas in doped semiconductors may undergo 
a ferromagnetic transition. The phenomenon occurs either at very low densities, due to the 
Coulomb interaction JtJ or, for instance, in GaAs under heavy doping with magnetic Mn 
impurities. As the ferromagnetic transition temperature T c is approached the longitudinal 
spin stiffness vanishes and so does the spin diffusion constant, which thus exhibits a critical 
behavior. 

We begin our analysis by assuming, as usual, a linear relationship between the number 
current densities J a (r) (a =t or J,) and the gradient of the local electro-chemical potentials 
ipa(r) = <t>(f) — iX I e )d f {n-\ 1 rij, T)/dn a , where cf>(r) is the electrostatic potential, e is the 
absolute value of the electron charge, and , ni, T) is the free energy per unit volume of a 
homogeneous interacting electron gas at the local spin densities n a (r) and temperature T ]To[| . 
This leads to the equation 

eJ a {r) = (a Q/3 (r)V0(r) - eD a p{r)Vnp{r)) (1) 

where <j a p is the homogeneous conductivity matrix, calculable from the Kubo formula and 
the diffusion matrix D a p is given by the generalized Einstein relation 

e 2 D a p = Y^Ua-fS-yp (2) 

7 

where 

^ _ d 2 f(n hni ,T) 
a>3 dn a dn/3 

is the static spin-stiffness matrix - minus the inverse of the spin susceptibility matrix. 

On a formal level the main effect of the Coulomb interaction is the appearance of non- 
vanishing off-diagonal elements of the conductivity and spin-stiffness matrices. These off- 
diagonal matrix elements have a simple physical interpretation, an ^ implies that an 
electric field acting only on the up-spin electrons must necessarily drag along a current 
of down-spin electrons. Similarly, S-^ ^ means that the chemical potential of up-spins 
df(n-\-,ni,T)/dni; is a function of both up and down spin densities. In addition the Coulomb 
interaction significantly modifies the values of the diagonal elements of these matrices, as we 
shall see momentarily. 

Let us apply eq. ([!]) to the calculation of the time evolution of a spin packet obtained 
by injecting an excess spin density Am(r, 0) = AMS(r) near the origin at time t = 0. We 
denote by m(f, t) — ?i|(r, i) — 7^x(r, t) the net spin density at point r and time t, by to*- ) = 
— the uniform value of the spin density at thermodynamic equilibrium, and by 
Am(f, t) = m(f, t) — m^ ' the excess spin density following spin injection. Similarly we denote 
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by n(r,t) = n^(f,t) + n^(r, t) the total electron density and by and by An(r, t) the 
equilibrium and the excess density respectively. 

The solution of the problem requires as additional inputs the continuity equations for the 
number and spin densities 

aA "« t > . -vsm ( 4) 



dt 

dAm(f,t) Am(f,t) 
dt t. 



V-J m (f,t), (5) 



where J m = J-t — is the spin current density and t s is the spin relaxation time, which is 
very long |,[|. 

In practice, following the procedure familiar in the theory of bipolar carrier packets pL 
we first combine eqs. (^), (||) and (Q) to eliminate the V • E term related to the Poisson 
equation, and then impose the local charge neutrality constraint |l2j 

An^r,t)=-A ni (f,t). (6) 

This yields the result 

dAm(r,t) Am(r,t) _ . ,_ . 

at t s 

+Li s E ■ VAm(r, t) (7) 

where 

(n T +n i )g T g x 

and 

a \ + <7 i 

are the effective mobility and diffusion constants , and E is an externally applied electric 
field. Eqs. (||) and (||) reduce to the expressions presented in |J in the noninteracting case. 
The constants a a and D a are presently given by 



D a = D aa -D a& . (10) 

The fact that the conductivities enter a spin symmetric combination while the diffusion 
constants are in a spin antisymmetric combination reflects the fact that the electrostatic field 
has the same sign for both spin components, while the density gradients have opposite signs 
(see®). 

The solution of eq. (|?]) in a homogeneous and isotropic liquid is 

Thus, we see that a Haynes-Shockley-type experiment can in principle determine Ms a- n d 
D s independently, provided that r s is sufficiently long. 

Let us now proceed to the calculation of Ms an d D s . The necessary inputs are the con- 
ductivity and the spin stiffness matrices. The conductivity matrix is best computed from 
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its inverse, namely the resistivity matrix p a p whose general structure is determined by the 
principle of Galilean invariance and Newton's third law. The explicit form of p a p, extracted 
from eq. (3) of Ref.( §]), is 

p=(^~^ PU m . PU n ). (12) 

Here r a are the combined electron-impurity and electron-phonon scattering times, usually of 
the order of 10~ 3 — 10~ 4 ns, m* is the effective mass of the carriers, and pn is the spin drag 
transresistivity calculated in Ref. ( ||). After computing a a p = [p^ 1 ]^ the diffusion matrix 
is straightforwardly obtained from eq. (||), and p s , D s are calculated from eqs. (g) and @. 
The algebra is greatly simplified by the reasonable assumption that the scattering times for 
the two spin components are not too different, i.e., t\ = rj. = td [ jljj . Under this condition, 
the result is 

^ = — ~ (!3) 



and 



D s = ^^—^, (14) 
e bcl-pn/PD 



where S = d 2 f(n, m, T)/dm 2 is the spin stiffness, S c = huTu/ An^ny is the Curie spin stiffness 
of an ideal classical gas, and pu = m* /ne 2 To is the Drude resistivity. 

Eq. ( |l3| ) tells us that the mobility of the packet is not explicitly modified by the Coulomb 
interaction and in fact coincides with the ordinary homogeneous mobility. Strictly speaking 
this result is only valid under the assumption that up-spin and down-spin electrons have equal 
mobilities and thus drift at the same speed in an applied electric field. Coulomb interactions, 
being Galilean- invariant, cannot change the total momentum of such a uniformly drifting 
electron gas. 

The situation is completely different for the diffusion constant. As the spin packet spreads 
out the up- and down-spin currents are directed in opposite directions and friction arises: for 
this reason the expression for D s contains the spin-drag resistivity as a factor that reduces the 
numerical value of D s . In addition, the Coulomb interaction together with the Pauli exclusion 
principle reduces the energy cost of spin-density fluctuations (i.e., the spin stiffness) further 
decreasing the rate of diffusion of a spin packet. 

Fig.l presents the necessary ingredients to calculate D s . Fig. la shows that pfi - a negative 
number - vanishes at low temperature as (T/Tp) 2 peaking at a temperature of the order 
of the Fermi temperature T F . As the inset illustrates, the prefactor 1/(1 — Ph/pd) (see 
eq. (|l4|)) displays a marked dependence on the sample mobility, increasing with the latter. 
Fig. lb shows S rescaled by its non- interacting value S n i and its behavior at the onset of the 
ferromagnetic instability. We evaluated S numerically starting from the work of Tanaka and 
Ichimaru jl5) where the free energy density is calculated as a function of temperature, density, 
and spin polarization, by means of a self-consistent integral equation approach that satisfies 
the thermodynamic sum rules. 

In fig. 2 we plot D s /D c , where D c = p s kBT/e is the classical noninteracting diffusion 
constant, for n-doped GaAs in a range of densities that are relevant to the experiments of 
Ref. The solid line corresponds to our full-interacting calculation, while the dashed line 
to the non-interacting case. We see that the interaction correction is quite significant, and 
reduces the value of D s as expected. Despite this reduction, D s is still considerably larger 
than D c consistent with experimental observations. These results show that the effect of the 
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reduced spin stiffness dominates at iow and intermediate temperature, while the spin drag 
contribution dominates at high temperature (see inset of fig. |^) . 

In the non-degenerate limit T >> Tp(n) we find that D s /D c approaches 1. In the non- 
interacting theory this limit is approached from above because the leading correction to 
the spin stiffness coming from the quantum kinetic energy is positive: S n i — S c [l + A^n/2\/2] 
where At is the De Broglie thermal wavelength at temperature T. In the interacting theory 
instead, the leading correction to D s is due to the spin Coulomb drag term and it is negative. 
In fact, for T » Z>(ra), Pu/pd ~ [n/(kBT) 3 / 2 ] ln(n/(fcsT) 2 ), and the logarithmic term 
dominates over corrections entering the spin stiffness in both the n — > and T — ► oo limits. 
Due to interactions, D s /D c — > 1 from below always. 

A very interesting feature of eq. (jlj) is the possibility of large variations in D s when the 
electron gas undergoes a ferromagnetic transition. From the curves in fig. lb and fig. 3 we see 
that S (interpreted as longitudinal spin stiffness in the ferromagnetic phase) and D s vanish at 
the transition temperature T c . For T < T c , D s /D c increases at first, due to the sharp increase 
in spin stiffness (sec fig.l), but then begins to saturate and tends to 1 as full polarization sets 
in II- 

In an ordinary electron liquid, the ferromagnetic transition is predicted to occur only at 
extremely low densities 0. There is, however, an interesting variant: semiconductors doped 
with magnetic impurities (Mn) can undergo a ferromagnetic transition at rather high carrier 
densities [p[|l7[|, n ~ 10 20 cm -3 for (Ga,Mn)As, and temperatures as large as 110K |J. Our 
theory can be extended to such materials, but, for the sake of clarity, we shall present this 
extension in a longer paper. Here we simply remark that in a simple mean field theory, such 
as that of Ref. ( |L8|| ), the calculated value of D s /D c for realistic values of the parameters 
has the form shown in fig. 3b, where, as in the intrinsic ferromagnetic case, D s vanishes at 
/• 0- 

In conclusion, we have demonstrated that a Haynes-Shockley experiment measuring D s 
and p s for a unipolar spin packet would be a sensitive probe of many-body effects such as 
the spin-Coulomb drag and the Coulomb enhancement of the spin susceptibility, and would 
provide a strong signature of a ferromagnetic ordering transition. Conversely, many-body 
effects must be taken into account in a quantitative theory of spin diffusion. 

* * * 
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Fig. 1 - (a) Transresistivity p-fi vs the reduced temperature T/Tf for typical semiconductor pa- 
rameters. The inset shows the behavior of the factor 1/(1 — Ph/pd) for three different mobilities: 
fi = 10 2 cm 2 /Vs (A), n = 3 x W A cm 2 /Vs (B), fj, = 10 4 cm 2 /Vs (C). (b) Spin stiffness S vs T/T F . The 
density is n — 4.2 x 10 11 cm -3 for the lower curve and increases by a factor 10 for each line starting 
from the bottom. The cusps represent the onset of ferromagnetism. 
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Fig. 2 - The interacting diffusion constant of a spin-packet (I) and its non-interacting approximation 
(NI) (rescaled by D c ) vs density for different temperatures. The inset shows the comparison with 
the value obtained considering interactions only through the spin Coulomb drag effect (D). In all 
the calculations the dielectric constant of the semiconductor is e = 12 and the mobility is fi — 
3 x 10 3 cm 2 /Vs. 
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Fig. 3 - (a) The diflFusion constant of a spin packet vs temperature in a low-density electron gas with 
no magnetic impurities, (b) Same as (a) for an electron gas of typical density n = 1.2 x 10 19 cm' 3 in 
a semiconductor doped with magnetic impurities (nj and Si are the impurity concentration and spin 
respectively). 



